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We report a theoretical analysis of the exclusive non-leptonic decays of the B and B° mesons 
into two light mesons, some of which have been measured recently by the CLEO collaboration. 
Our analysis is carried out in the context of an effective Hamiltonian based on the Standard 
Model (SM), using next-to- leading order perturbative QCD calculations. We explicitly take into 
account the 0(a s ) penguin- loop diagrams of all four- Fermi operators and the 0(a s ) tree- level 
diagram of the chromomagnetic dipole operator, and give a prescription for including their 
| effects in non-leptonic two-body decays. Using a factorization ansatz for the hadronic ma- 
trix elements, we show that existing data, in particular the branching ratios B(B ± — ► rj'K^), 
B(B ± -> ^K ), B(B°(W) -> and B{B ± -> wh*) (/i* = n ± ,K ± ), can be accounted 

for in this approach. Thus, theoretical scenarios with a substantially enhanced Wilson coef- 
ficient of the chromomagnetic dipole operator (as compared to the SM) and/or those with a 
substantial color-singlet cc component in the wave function of rf are not required by these data. 
We predict, among other decay rates, the branching ratios for the decays B°(B°) — > 7r ± 7r =F and 
B^ — > 7r°7r ± , which are close to the present experimental limits. Implications of some of these 
measurements for the parameters of the CKM matrix are presented. 
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1 Introduction 



Recently, the CLEO collaboration has reported first measurements in a number of exclusive 
decays, B ± -> r/K*, B±_^ ^K , B°(W) -> ir^K*, B ± -> ^/i* (/i* = ?r ± , i^ ± ), -> 
K°/i ± , £± -> vr /^, -> /^tt* and the inclusive decay 5± -»• 7/ + X || |, |,_which 

involve the so-called QCD penguins. In addition, a number of related decays such as B°(B°) — > 
7r T 7r ± and B ± — > 71-^71-° are on the verge of measurement 0]. On the theoretical side, considerable 
effort has gone into studies of non-leptonic weak decays in terms of estimating decay rates M 



- 13 and the inherent direct and indirect CP asymmetries | 1Q| , fT4| - [p2| . Since the first 
measurements of the above-mentioned decays, theoretical interest in this subject has surged 
and recent literature is rife with all kinds of interesting interpretations of data, both within 
and beyond the SM [2vJ - |^9| . Of these, the decay mode B ± — ► r/K ± is conspicuous due to its 
reported high branching ratio B{B ± -> r/K*) = (7.1±|;f ± 0.9) x 10~ 5 [g[. 

The standard theoretical framework to study non-leptonic B decays is based on the effective 
Hamiltonian approach, which allows to separate the short- and long-distance contributions in 
these decays using the Wilson operator product expansion [pOR . QCD perturbation theory is 
then used in deriving the renormalization group improved short- distance contributions [31 



This program has now been carried out up to and including the next-to-leading order terms 



32j p3| , but the long-distance part in the two-body hadronic decays B — > M\M 2 involves the 
transition matrix elements {MxM 2 \Oi\B) at a typical hadronic scale, where Oi is an operator in 
the effective Hamiltonian (see below). Calculating these matrix elements from first principle is 
a true challenge in theory which remains to be met. In view of this, a number of approximate 



schemes has been put forward. The one we use here is based on the idea of factorization [34, 35 



in which the final state interactions (FSI) are assumed absent, and hence the hadronic matrix 
elements in the decay B — > M\M 2 factorize into a product of two comparatively more tractable 
matrix elements. These are then taken either from data or calculated in well defined theoretical 
contexts, such as QCD sum rules and potential models jjj |13j, |36|]-[[4(|. This framework does 
remarkably well in accounting for non-leptonic two-body B decays involving the current-current 
operators 0\ 2 @, 0] (see section 2 for definition). Recent analyses have shown that data on two- 
body non-leptonic B decays on the so-called heavy-to-heavy transitions, such as B — > (D, D*)h, 
B — ► (D s , D*)D, B — ► J/iph (with h being a light hadron), can be described in terms of 



two phenomenological parameters, a\ and a 2 m, whose values seem to be universal [41, 13 



Techniques based on HQET [H2J allow in some limited cases to "derive" such factorization 



properties [^3| and yield results which are in agreement with data. 

Motivated by the phenomenological success of factorization in the heavy-to-heavy non- 
leptonic B decays, we would like to pursue further this framework in the domain of the so- 
called heavy-to- light transitions, B — > h 1 h 2 , where hi and h 2 are light hadrons. The recently 
measured B decays 0, [| belong to this category and they should be analyzed on their own, 
without prejudice about the suggested values of the effective parameters from the heavy-to- 
heavy transitions. The decays B — > h\h 2 in most cases involve mixing among the current- 
current, QCD-penguin and the chromomagnetic operators. Our hope is that once perturbative 
QCD corrections are taken into account, these decays may allow themselves to be described in 
terms of a few phenomenological parameters. Related work along these lines concerning QCD 
penguins in non-leptonic B decays has been done prior to this analysis || 10], which we make 
use of here, improve upon and extend. 
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Our analysis is based on the following three main ingredients: 



We work at next-to- leading logarithmic (NLL) precision, taking into account the 0(a s ) 
one-loop penguin-type diagrams of all four-Fermi operators in the effective Hamiltonian 
and some process-independent parts of the vertex correction diagrams associated with 
these four- Fermi operators. We also take into account the effect of the 0(a s ) tree-level 
matrix element associated with the chromomagnetic dipole operator via the process b — > 
sg -> sq'q'. 

To calculate the hadronic matrix elements, we propose a simple factorization ansatz which 
allows to include the effects of the 0(a s ) matrix elements just discussed above. 

In calculating B decays involving an 7/ or 77 meson, such as B — > rj'K ± and B — > r]K ± ) 



we include the contribution from the decay b — > s(cc) — > s{j],rf) ||4], (45| . The required 



decay constants and mixing parameters are estimated using data on the radiative decays 



J/ip — > r/ c 7, 7/7, 777 and the two-photon decays of the 77, 77', and r] c |j46[ . Concerning (77, 77')- 
mixing, we discuss both the conventional (one mixing-angle) formalism [[Tj], ^8] and the 
one involving two mixing-angles in this sector, which is suggested by the 1/iVc-improved 
U(3) ® U(3) chiral perturbation theory framework ^{J 50]. Since the latter formalism is 



also favoured by a recent phenomenological analysis |51] of the data on the 777 and 77^ 



form factors [|52|-||55||, we use it in our estimates for B decays involving 77' and 77-meson. 

We would like to make a number of remarks pointing out the overlaps and differences with 
earlier analyses and explaining our factorization ansatz. 

Concerning the QCD-perturbative part, we note that our calculations come close to the 
derivation given in |10[ but are more complete as far as the NLL contribution is concerned. 



We find that the NLL improvements implemented by us reduce the scale dependence in var- 
ious non-leptonic decay rates. This result is in line with what has been demonstrated in the 
radiative decays B —>■ X s + 7 in the same accuracy . Further, the complete NLL contribu- 
tion is important numerically, both compared to the leading order result and the NLL result 
obtained by keeping only the charm penguin contributions from the operators 0\ 2 - We show 
this quantitatively in the context of the branching ratio B(B ± —>■ Kir ), comparing it with the 



estimates of the same based on keeping only the 0\ 2 penguins [^] . 

Concerning the second point noted above, we remark that our factorization prescription 
introduces just one free parameter, called £, which is supposed to compensate for the neglect 
of color octet-octet contribution in evaluating the hadronic matrix elements in the heavy-to- 
light sector B — > h 1 h 2 . This modifies the strength of the effective coefficients ai,...,as from 
their perturbatively calculated values (see section 3). Clearly, this is the simplest ansatz and 
may have to be modified eventually as more precise data on heavy to light B decays become 
available. 

We discuss the last point mentioned above concerning the decays B ± -> rf{K ± ,K* ± ) and 
B ± —>■ r](K ± ,K* ± ). Expressing the charm quark content in the 77' meson in terms of the 
matrix element (77' 1 07^750(0) = —if^ q^ we find using data on the J/tp 7/ c 7 and J/ip — > 77^ 
decays that |/i | — 5.8 MeV. The corresponding decay constant for 77 meson is estimated to 
be |/^ c '| — 2.3 MeV in the conventional (77, 77')-mixing formalism and \fjj C '\ — 0.93 MeV in the 
l/iV c -improved approach. The decrease in the value of \f^\ reflects the small value of the 
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singlet mixing angle 8q, which makes the rj an almost pure octet state [|{J, hence also reducing 
the cc component of the 77-meson. Our estimate |/~? | — 5.8 MeV is to be contrasted with the 



range = (50 — 180) MeV obtained in |23[. (Likewise, we find \f!% / frf\ — 0.08, which is also 
an order of magnitude smaller than the one given in pi}] .) We note that our estimate of \f!^\ 
is consistent with the bounds —65 MeV < < 15 MeV, which have been obtained in the 
meanwhile from an analysis of the Q 2 -dependence of the electromagnetic form factor of rf |5l 



Likewise, data on the electromagnetic form factor of 77 is consistent with |/^| being small [51 
With our estimate of |/,^|, we find that this charm-induced contribution does not dominate 
the matrix element for B ± — > rj'K ± ; the penguins play a more important role numerically in 
this decay. 

The branching ratio B{B ± — > rfK ) as well as those of the related ones B — > rfK , 
B ± — > t]K ± , and B^ — > rjK^, depend upon the interference of the amplitudes arising from the 
chain b — > s{cc) — > s(f]',r]), and the ones arising from calculating the matrix elements of the 
rest of the operators. Concentrating on the decay B — > rfK , we note that the sign of the 
term involving the b —>■ s(cc) — > s(r]', rj) in the full amplitude is not determined a priori. Since 
the solutions with constructive or destructive interference terms are both logical possibilities, 
we have estimated B(B ± — > rfK^) for both cases, with the positive-/^ solution yielding 
marginally larger rate. However, more importantly, we find that the rate in this decay (and in 
some others) depends significantly on the parameter £. Hence, to make absolute predictions, 
the phenomenological value of this parameter has to be determined. We study a number of 
measured B — > h\h2 decays to estimate a range for £ which, given the present experimental 
errors and theoretical accuracy of our approach estimated by us as a factor 2 in rates, is 
understandably not very precise at this stage. The range < £ < 0.5 is consistent with data. 

This paper is organized as follows. In section 2, we review the effective Hamiltonian for the 
non-leptonic B decays and calculate the matrix elements of the operators at the quark level in 
the NLL precision. In section 3, we formulate our factorization ansatz to calculate the hadronic 
matrix elements in the two-body decays B — > h\h2- The matrix elements for various decay 
modes of interest are also detailed here, together with a brief review of the mixing formalism 
for the r\ - rf — r) c sector. Our estimates of the decay constants f& and ffl relevant for the 
decays B ± -> (rj 1 ,rj)(K ± , K #± ) are also given here. Section 4 contains our numerical results. 
The input values for the various quantities (coupling constants, form factors, quark masses) 
are collected here in several tables. We compare the branching ratios with the CLEO data 



varying the factorization-related parameter £ and parameters of the CKM matrix [|57 |. The 



potential impact of some of these decays on the CKM phenomenology is illustrated in terms of 
the ratios of the branching ratios, which are more reliably calculable. In particular, the ratios 
Rx = B(B (W) -> n^K^/BiB^ -> ^K) and R 2 = B(B°(W) -> ^h^/B^ -> tt^K), 
constrain the CKM-Wolfenstein parameter p and rj ||58|| . The potential importance of R± in 
determining the angle 7 has been emphasized by Fleischer and Mannel |2P[. Interestingly, 



within the theoretical framework presented here, the measured ratio R\ = 0.65 ± 0.40 suggests 
(at ±lcx) that p > 0, which in turn implies 7 < 90°, where 7 is one of the CP- violating angles 
of the unitarity triangle. We also comment on the effect of an (assumed) enhanced coefficient 
of the chromomagnetic operator, Cs(mw), in non-leptonic two-body B decays. This scenario 
has been discussed in the context of new physics effects in B decays |59], |60|, . We find, using 
the decay B ± — > K^, that varying the ratio C%{mw) /C%(rriw) SM in a large range (±10) has 
no appreciable effect on the branching ratio within the present accuracy. Finally, we conclude 
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with a summary in section 5. 



2 Effective Hamiltonian for the non-leptonic decays B 

hih 2 

We write the effective Hamiltonian i? e // for the AB = 1 transitions as 



H eff = — 



Gp 

71 



v ub v: Q (do? + c 2 o u 2 ) + v cb v; q {c x o\ + c 2 o c 2 ) - v tb v* £ao 



i=3 



(2.1) 



where q = d,s and Cj are the Wilson coefficients evaluated at the renormalization scale /i; the 
current-current operators 0"' c and G^ c read 



2 " 



{u a b a ) v _ A ( y q f 3Uf 3 ) v _ A 



(uf3b a ) v _ A (q a Up) v _ A 
while the QCD penguin operators O3 — Oq are 
3 = (q a ba) V -A E 

Finally, the dipole operators O7 and 0% read 



9 c 



o 5 

Oa 



(cf3b a ) v _ A (q a Cf3) v _ A 



(2.2) 



{qpb 



a)V-A 



(2.3) 



7 
8 



(e/87r 2 )m 6 sa^ (1 + 75) 6 F. 



(^/8tt 2 ) m 6 s a (1 + 75) (A^/2) G 



(2.4) 



Here a and /3 are the £77(3) color indices and A^g,A = 1,...,8, are the Gell-Mann matrices. 
The subscripts V ± A represent the chiral projections 1 ± 75. Thus, in eqs. (^.2[) and 



(iXafysV-.A = m q 7 m (1 — 7s) fys etc. In eq. (|2.4|) and G^ u denote the photonic and the gluonic 
field strength tensor, respectively. We note that we neglect the effects of the operator O-j in 
the present analysis as well as the so-called electro-weak penguin (4-Fermi) operators which 
we did not list explicitly. Likewise, the effect of weak annihilation and exchange diagrams will 
be neglected here. This is in line with the investigations reported in literature [[HJ. Working 
consistently to NLL precision, the coefficients Ci-Cq are needed in NLL precision, while it is 
sufficient to use the LL value for C§. The relevant Wilson coefficients to the desired accuracy are 
listed in table 1 for the two scales fi = 5.0 GeV and \i = 2.5 GeV, where Cy = C 7 — C 5 /3 — C 6 
and Ct ff = C 8 + C 5 . 



2.1 Quark- level matrix elements 

The Wilson coefficients of the four-Fermi operators depend on the renormalization scale; in 
addition, in NLL precision, they also depend on the renormalization scheme. These unphysical 
dependences are compensated in principle by a corresponding scheme/scale dependence of the 
matrix elements of the operators. However, when using the factorization ansatz for the hadronic 
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fj, = 5.0 GeV 


H = 2.5 GeV 


qNLL 


1.070 


1.117 


qNLL 


-0.166 


-0.257 




a ni i 

0.011 


0.017 


(jNLL 


-0.031 


-0.044 


qNNL 


0.009 


0.011 




-0.037 


-0.056 


(jeff,LL 


-0.303 


-0.338 


r eff,LL 


-0.144 


-0.158 



Table 1: Wilson coefficients CA^i) at the renormalization scale /x = 5.0 GeV and /j = 2.5 GeV 
in the NDR scheme, d -C 6 are in NLL accuracy, while C e 7 ff and C 8 e// are in LL precision. For 
a s (/i) (in the MS scheme) we used the two-loop expression with 5 flavors and a^ s {mz) = 0.118; 
mf s (m t ) = 165 GeV (equivalent to m tiPO i e = 175 GeV).. 

matrix elements of the operators, these cancellations do not take place in practice, because the 
factorized matrix elements of the operators are expressed in terms of decay constants and form 
factors, and are as such scheme/scale independent. To achieve this cancellation, we include 
perturbative QCD corrections to the partonic matrix element before doing the factorization 
step. We fully calculate the one-loop penguin-like diagrams in Fig. [l|(a) and some process 
independent parts (see below) of the vertex correction diagrams associated with the four- Fermi 
operators, as shown in Figs. 0(a). These two classes of corrections are sufficient concerning the 




Figure 1: (a) Penguin-type diagrams associated with the four- Fermi operators Oj (i = 1, 6); 
(b) Tree-level matrix element of the chromomagnetic dipole operator Og. 

cancellation of the scheme/scale dependences. Furthermore, the contribution associated with 
the operator 8 , where the gluon splits into a quark - antiquark pair, as shown in Fig. |l|(b) is of 
the same order in a s as the corrections just mentioned and is therefore also taken into account 
in our analysis. 

As we use in this paper the Wilson coefficients obtained in the naive dimensional scheme 
(NDR) with ant i- commuting 75, we also have to evaluate the various 0(a s ) corrections in this 
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(a) 



(b) 



Figure 2: (a) Vertex correction diagrams to the four-Fermi operators O, (i = 1, ...,6); 
(b) Corresponding Bremsstrahlung corrections. 



scheme. These corrections can be absorbed into effective Wilson coefficients C^, which for a 
general SU(N) color group can be written as 



eff 



C{ 



c% fs 



c e J f 



Ci + 

c 2 + 
c 3 - 



Air 

Air 
1 



ry + ly lo S 



C,+ 



+ 7y log — I C,- + 

V ^ / 2j 



a u 



2NA 



7T 



(C t + c p 



C4 + ^(a + C p + C 9 ) + ^ ^ + 7 ^log^ c j + 



2A^4vr 



(C t + C p + C s ) + -^ 



-[r v + lv \o S -^C 1 + 



Ce + ^ (C( + Cp + Ca) + ^( r j + 7 y log ^Q + 



(2.5) 



We have separated the contributions Ct, C p , and C g arising from the penguin-type diagrams of 
the current-current operators 0i,2, the penguin-type diagrams of the operators O3-O6, and the 
tree- level diagram of the dipole operator Og, respectively. The process-independent contribu- 
tions from the vertex-type diagrams are contained in the matrices ry and jy. Here jy is that 
part of the anomalous matrix which is due to the vertex (and self-energy) corrections. This 
part can be easily extracted from 7^ in ref. [33[|: 



lv 



( "2 
6 




V 



6 

-2 










-2 
6 







6 

-2 









2 




\ 




-6 
-16 ) 



(2.6) 



The matrix ry contains constant, i.e., momentum-independent parts associated with the vertex 
diagrams. This matrix can be extracted from the matrix f defined in eq. (2.12) (and given 
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explicitly in eq. (4.6) in ref. |33[] ) : 



rv 



( 


7 
3 


-7 














\ 




-7 


7 
3 
























63 
27 


63 
7 9 


















-7 


3 
























1 

3 


1 




V 














-3 


35 
3 


/ 



(2.7) 



Note that the /z dependence and the scheme dependence of the vertex correction diagrams are 
fully taken into account in eq. (|2.5| ) by the terms involving the matrices 7y and ry, respectively. 
There are, however, still scheme-independent, process-specific terms omitted as indicated by 
the ellipses. When calculating inclusive quantities, such as the semileptonic branching ratios 
and .B-hadron lifetimes, it is straightforward how to take these corrections into account. The 
virtual corrections are infrared divergent on their own, but together with the Bremsstrahlung 
contributions in Fig. |2](b), they lead to a finite and well defined 0(a s ) correction, which is found 
to be small. However, it is less obvious how to include them in exclusive two-body decays. The 
point is that the division of the final states with and without the extra gluon is ambiguous and 
can be meaningfully defined only with a cut-off. As such a separation into virtual corrections 
and soft gluon Bremsstrahlung contributions is arbitrary anyhow, we only take into account the 
terms involving 7y and ry. The explicit 0(a s ) contributions which emerge from the penguin 
operators involving a qq pair in the loop are infrared finite on their own and hence do not 
require a cut-off. 

The quantities C t , C p , and C g given by the diagrams shown in Fig. [1] read in the NDR 
scheme (after MS renormalization) 



C t = -Ct 



V V* V V* 



V tb V t * q 



v tb v t * q 



CAm) 



2 
3 



2 , m 2 
3 l0 V 



?7r 



(2.8) 



Cn — C-t 




+(c A + c 6 ) E 

i=u,d,s,c,b 

c n = 



2 , mj A „ / a 

3 '^-^l- 



2m h 



(I 2 ) 



with C 8 e// 



mt 



C*8 + C*5- The function AFi(z) is defined as 

AFi(z) = -4 / dxx(l-x) log[l-zx(l 
Jo 



ie\ 



(2.9) 



(2.10) 



(2.11) 



Two remarks are in order here. First, our expressions for G\ in eq. (|2T5| ) are written in 
terms of the Wilson coefficients in the NDR scheme. Analogous expressions (but with ry = 0, 
and C n = 0) have been obtained earlier in the literature [Hj. Comparing the expressions 



Iv 



given here with the ones in fllQ |, where the corresponding quantities c% , c e J^ are expressed 
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in terms of the so-called renormalization-scheme-independent Wilson coefficients q introduced 
in ref. |33], one notices that the constant terms appearing explicitly in C t and C p in the two 



papers are different. As the scheme dependence cancels automatically when including the one- 
loop matrix elements discussed above, we prefer to work with the Wilson coefficients in the 
NDR scheme. 

Second, we have to explain the assumption which allows us to absorb the tree-level diagram 
b — > sg —>■ sq'q' associated with the operator 8 into the contribution C g appearing in the 
expressions for . It is straightforward to write down the matrix element 

(sq'q'\Os\b) = -^^ (^(1 + 75)^) , (2-12) 

where q is the momentum transferred by the gluon to the (g',g')-pair. In the factorization model 
to be described below, q' and q' cannot go into the same meson in the process B — > h\hi due 
to color, i.e., q' goes into hi, while q' goes into h 2 or vice versa (see Fig. |](b)). The quantities 
Ct, C p and C g depend on the momentum q. Since we are interested here only in two-body 
decays, we assume for simplicity, that the three-momenta of q' and q' are equal in magnitude 
but opposite in direction in the rest frame of the 6-quark. The momentum transfer q is then 
proportional to pb, i.e., 

9" = ^^ , (2-13) 

1 1 Of) 

where (q 2 ) is an averaged value of q 2 . Inserting ( |2.13| ) into eq. ( |2.12| ) and using the equations 
of motion, the expression for C g in eq. ( |2.10|) is readily obtained. To be consistent, we should 



also replace q 2 by (q 2 ) in the expressions for C t and C p in eqs. fl2.8|) and (|2.9|) , respectively. 
To estimate the theoretical uncertainty introduced thereby, we treat (q 2 ) as a parameter which 



varies in the range < (q 2 ) < ml/2, following the prescriptions in literature \ Wj , |I~7 



To summarize: The various 0(a s ) corrections have been absorbed into effective Wilson 
coefficients C\ * {% = 1, 6); these coefficients are scheme independent and the term ~ a s log \x, 
which dominates the scale dependence of the original Wilson coefficients C, and the one-loop 
matrix elements, is absent in C±. What remains to be done is to estimate the hadronic matrix 
elements (Ji\h-i\C^ Oi\B) for i — 1, 6. The numerical values of the quantities are given 
in table 5 in section 4.1.4. 



3 Factorization ansatz for the matrix elements in B —* 

hih 2 

We have now to work out the hadronic matrix elements of the operators Oj (z = 1, 6) for the 
processes of interest. We use the factorization approximation, which we briefly explain for a 
specific example. Consider the matrix element due to the u-quark contribution of the operator 
O5 for the process B~ — > K~u, i.e., 

(K-u\0^\B-) , Of = (*7m (1 - 75)6) {ui» (1 + ls)u) . (3.1) 
There are two contributing diagrams D\ and D2 shown in Fig. |^. The factorization approxi- 
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CO 




K" 




D, 



Figure 3: D\ and D 2 are the two diagrams contributing in the factorization approximation. 
See text. 



mation for D\ is readily obtained: 

D x = (u\u-f (1 + 75) w|0> (^-1^(1-75)615-) = {u\uu-\0) (K-\sb„\B~) , (3.2) 

where here and in the following the short-hand notation qq'__ stands for 

QQ- = tin (1 - 7s) q ■ (3.3) 

To get D2 in the factorization approximation, we first write the operator in its Fierzed 
form 



O 



(«) 



-2 (up(l - 7 5 )6 Q ) (s a (l + 75)14/3) 
= -2 



i («(1 - 75)6) (5(1 + %)u) + \ (u{l - 75 )A6) (s(l + 75 )Am) 



where A denotes a color matrix. Only the first term in the square bracket in eq. 
color singlet- singlet) contributes in the factorization approximation. One gets 



D 2 = - (u\u(l - j5)b\B~) (K-\s(l + 75 )«|0) 



using the Dirac equation, we can write D 2 as 



Do = - — 



K 



{u\ub_\B-) (K-\su_\0) . 



(3.4) 
(being 

(3.5) 
(3.6) 



N (m s + m u ) (m b + m u 
Doing analogous manipulations, the complete matrix element M for B~ — > uK~ , defined as 



M = (uK-\H eff \B-) ; H eff = ^= 



V U bV* s (ct ff o\ + c e 2 ff oi) - V tb V* j2 Cf f Oi 



i=3 



(3.7) 
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is then easily obtained. One gets 

(ai (K~\su-\0) (u\vb-\B~) + a 2 (K~\sb-\B~) (uj\uu-\0) 
2a e m 2 K 



M 



-v*K 



a 4 



(K-\su_\Q) {u\ub_\B-) 



The quantities a 
in eq. (gj): 



+2(a 3 + a 5 )(X~|s6-|5-}(a;|«w_|0))} . (3.8) 
6) are the following combinations of the effective Wilson coefficients 



&2i-l 



2i-l T TT ) a 2i — <-^2i + 



°2i-l 



1,2,3 



^y- Z2 7 ZZ 1 yy ZZ — 1 7 7 7 

(3.9) 

The explicit 1/N terms in eq. ( |3.9| ) are always accompanied by an octet-octet contribution; 
this can be seen explicitly in eq. (|3.4j ). As one discards this octet-octet contribution in the 
factorization approximation, one usually replaces 1/N by £ and treats £ as a free parameter 
with the hope to compensate phenomenologically for the omitted octet-octet contribution in 
terms of a rescaled value of £. Note, however, that the 1/N factors appearing explicitly in the 
(perturbative) expressions for the effective Wilson coefficients in eq. ( |2.5|) are not replaced by £ 
in our work in contrast to |]IU|| , where also these 1/N factors were replaced by £. We think that 
replacing also these 1/N terms by £ destroys the scheme independence of the effective Wilson 
coefficients. 

It is worth pointing out that the factorization Ansatz just discussed is the simplest one. 
Also, it is implicitly assumed that the relative strong phases (such as the ones arising from 
the final state interactions but also due to the non-perturbative contributions to the charm- 
penguins), contributing to the different a^s are small. Of course, it does not mean that the 
strong interaction phases are assumed absent. The ones generated by the next-to- leading order 
perturbative QCD contributions from the charm penguins are taken into account. It remains 
to be seen if the non-perturbative phases from the competing tree and penguin contributions 
in the processes discussed here are indeed small. 

Finally, before giving the matrix elements for the various exclusive two-body decays, we dis- 
cuss the parametrization of the decay constants and form factors which appear in the factorized 
form of the hadronic matrix elements. The form factors are parametrized as 



{P(p')\V,\B(p)) = [(p + p') 



m B — m\ 



2 2 

m B — nip 



q»F (q 2 



(3.10) 



(V(e,p')\(V,-A,)\B(p)) = 

VTiy 



—{m B + m v 

e* ■ q 
m B + m v 



e* ■ q 



A 1 (q 2 ) + 



(P + P) 



m B -m v 



A 2 (q 2 



e -q—z 



q^Aoiq 2 



where P(V) is a pseudoscalar (vector) meson, q = p — p', 



MO) = ^^Mo) -^^A 2 (o) 



2m, 



2m 



v 



(3.11) 



(3.12) 
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and -Fi(O) = F o (0). The decay constants f P and f v are defined as 

(0\A li \P(p))=if PPli , (0\V li \V(e,p))=if v m v e li . (3.13) 

With these definitions we are in a position to write down the formulas for the matrix elements 
for the two body decays. They are given below explicitly for the four generic decay modes: 
B — > 7T7T, B — > Kit, B^ — > K ± u,tv ± u and B ± — > (K ± , K^)^, rf), which are also the ones 
we calculate numerically in the next section. However, the formalism given here is general and 
applicable to all two-body S-decays of the type B — > PP, B — > PV, and B — > VV. 



3.1 B -> 7T7T 

In this section we discuss the processes -B 1 * 1 — > vr 1 ^ , B°(B°) — > 7r ± 7r T and B°(B°) — > 7r°7r°. 



3.1.1 5 ± 



7T 7T 



The matrix element M for S 
SU(2) breaking effects) 



7r 7T° involves the operators 0" and O2 an d reads (neglecting 



1/ -j=V ub V: d {a 1 + a 2 ) (7r-|du_|0) <7r°|^6_|^- 



with 



(n-\du_\0) (7T°\ub_\B-) = iU (m| - m%) F W '«) 
The branching ratio B(B~ — > -K~n°) is then given by the expression 



PI 



B(B~ -> 7T-7r u ) = r B — |M| 2 



8tt 



(3.14) 



(3.15) 



(3.16) 



where tb is the lifetime of the 5°-meson and \p\ is the absolute value of the 3-momentum 
of the 7r~ (or the 7r°) in the rest frame of the B° meson. This expression for the branching 
ratio holds for other two-body decays being discussed with obvious changes of the indicated 
quantities. Hence, we shall give subsequently only the matrix elements M. Also, we shall give 
only the amplitudes for the decays of B~ and B°, and the matrix elements for the charge- 
conjugate processes are then obtained by complex conjugating the CKM factors. Since we are 
not addressing the question of CP-violation in this paper, all decay rates given later are to be 
interpreted in terms of the averaged branching ratios. Thus, for example, the branching ratio 
B(B ± — > tt^tt ) is defined as 



B(B 



± 



^71°) = \- (B{B- -> 71-71°) + B(B + -> 7T + 7T )) 



(3.17) 



3.1.2 B°(B°) -> tt+tt- 

The matrix element M for B° — > 7r~7r + reads 

V2 



M 



VubV*^ -v tb v; d U 4 + 



2a 6 ml 



(m h - m u ) (m u + m d ) 



with 



(7r _ |d«_|0) (n + \ub_\B°) , 
(7r-\du-\0) (7r + \ub_\B ) = if w (m 2 B 



m 



73 — >7T~ 



[m„ 



X 



(3.18) 
(3.19) 
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3.1.3 B°(B°) -> 

The matrix element M for 23° — > 7r°7r° reads 



M 



a 6 m£ 



m rf (m b - m d ) 



x 



with 



(7r |««_|0> (Adb-\B°) = ^ (m| - ml) F^\ml) . 



(3.20) 



(3.21) 



When calculating the decay width, we have to take into account an extra factor 1/2 due to the 
two identical particles in the final state. 

3.2 B — ► K7T) B — > ifiY 
3.2.1 B ± -> Ktt ± 



The matrix element M for .B 



7r K° reads 



M = 



V2 



v tb v* 



(I4 



2a 6 m 



K 



(m b - m d ) (m s + m d ) 



{n-\db-\B-) (K°\sd-\0) 



with 



(n-\db-\B~) (K°\7sd-\0) =if K (m 2 B - ml) F*^~ (m 2 K ) 



(3.22) 
(3.23) 



3.2.2 B°(B°) -> K^~7r~^ 

The matrix element M for B° — > tt+K - reads 



M 



KfeC «i " W t ! a 4 + 2a 6 



V2 

(Tr+lMfe^l^ ) (K~\su_\0) 



m 



K 



(m b - m u ) (m s + m u )) 



x 



with 



(vr + |ii 6_ | B°) {K-\su-\0)=i f K (m% - ml) (m 2 K ) 



(3.24) 
(3.25) 



3.2.3 B ± -> J^vr 

The matrix element M for — > K~tt is given by 



M 



v/2 



KbK!, ai + a 2 



m| - rr& F<f ^ K (ml) ' 



( m 
- VtbV ^ { a * + 2a \m b - mu )( ms + m u )) 



V2f K ml -ml F^°(m\)) 

(7r°\ub_\B-) (K-\su-\0) 



with 



{K-\su-\0) (K°\vb-\B-) =if K (ml-ml)F B ^(m 2 K ) , 

12 



(3.26) 
(3.27) 



3.2.4 B ± -> K°K± 

The matrix element M for — > K°K~ is given by 



M 



m 



A' 



(K-\sb-\B-)(K°\ds-\a) (3.28) 



with 



(Ky S .\0) (K-\sb-\B-) = if K (m 2 B - m 2 K ) F^ K {m 2 K ) . 



(3.29) 



3.3 B ± ^K ± u,tt ± u 
3.3.1 -> K±u 

The matrix element M for B~ — > reads 



ai + a 2 



B^K' 



(m 



A^(m 2 



v tb v* 



+a 4 - 



2qq m 2 K 



(m b + m u ) (m s + m u ) 



k) VZ/k. 
(u\ub-\B-) (K-\su-\0) 



2(a 3 + a 5 



(m 



f. 



A^(m 2 K ) V2f K 

(3.30) 



with 



(u\ub.\B-) (K-\su-\0) = -i/ Jf 2m w (p B .£*)^K) . 



(3.31) 



There is only one non-vanishing helicity amplitude. In the rest frame of the decaying B meson 
only longitudinally polarized u;'s are produced, ps • e* is then given by 



m B | | 

Pb ■ t w = \p\ , 



(3.32) 



where \p\ is the absolute value of the 3-momentum of the u (or the K ) in the B rest frame. 
3.3.2 B ± -> n^u 

The matrix element M for B~ — > h~oj reads 



+a 4 — 



ai + a 2 



2a 6 m£ 



(2a 3 + a 4 + 2a 5 ^ 



(W|U6_|S") (7T _ |du_|0) 



with 



(ra 6 + m u ) (m d + m u )_ 
(cu\ub_\B-) (7r-|d«_|0> = -i /,2m, {p B ■ el) <^K) 



^K) y/2f n 

(3.33) 
(3.34) 



3.4 Mixing in the r\—r(—r\ c system and the decays B ± — > K ± rj^'\ K* ±r r]^ 

Before we write the matrix elements for B ± — > K ± r]', B ± — > K* ±r q' ', B ± — > and 5 ± — > 

j^*±^ j n f ac torization approximation, we give a short discussion about the rj - rj' — rj c 
system. Our main interest for the decays mentioned above is to compute the hadronic matrix 
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elements (0 | (07^75) | 77) and (0| (07^75) \r]'). The conventional (77, ^-mixing formalism involves a 
single mixing angle (called henceforth 9) and it has been argued that it provides a satisfactory 
description of the decays involving rj and rj' [|7 , f48|| . However, recently the inadequacy of this 
mixing formalism has been pointed out in the context of the l/iV c -improved chiral U(3) <g> U(3) 
perturbation theory [^9], pD| . Instead, a formalism which involves two mixing angles in the 
5 , f/(3)-octet and -singlet sector (called henceforth 9 8 and #0) is proposed. Since, the SU(3)- 
singlet component \rj ) in general mixes with the |cc) component, introducing another angle 
(called 9 c5 ), we shall term the resulting mixing formalisms as the two-angle (involving 9 and 
9 c5 ) and three- angle (involving 9 , 9$ and 9 cS ) frameworks. 



3.4.1 rj - if — T) c system in the two-angle mixing formalism 

Here, the physical 77 and 7/ states are considered as mixtures of the rjs and rjo states (47 

|?7) = cos^lr/g) — sin^^o) , \rf) = sin 9\r]$) + cos 9\r] ) , 



(3.35) 



where r] 8 belongs to the SU(3) octet of the light pseudoscalar (Goldstone) mesons, while rj is 
an SU(3) singlet. In the quark basis they are given by 



\Vs) 



1 

7q 



\uu + dd — 2ss) , 1 770) 



1 

7s 



\uu + dd + ss) 



(3.36) 



The mixing angle 9 can be extracted from the measured ratios |47 



and 



Y(n -> 77) 

r(7r° — > 77) 



r(V -» 77) 
r(7r° — > 77) 



J 7T 



cos el + e 2 d - 2e 2 sin el + e 2 , + e 2 



fo 



(3.37) 



f 2 



sin 61 e 2 + e 2 , - 2e 2 cos # e 2 + ef, + el 



v/6 



+ 



/o 



v/3 



(3.38) 



where ej are the quark charges, / 8 an d fo are the decay constants of the pion, the eighth 
component of the octet, and the singlet, respectively. Using fs/ f w = 1-34 ± 0.03 |J2] and the 
measured decay widths @ r(7r° -> 77) = (7.7 ± 0.55) eV, T(r/ -> 77) = (0.46 ± 0.04) KeV, 
W -> 77) = (4.26 ± 0.19) KeV, one obtains 



= -21.3° ±2.5° 



1.09 ±0.05 



(3.39) 



It has been suggested in the context of the radiative decays J/i/j —> 777, 7/7 that they can be 
enacted by modelling them on the decay chain J/ip — > ^7 — > (77, t)')j, involving the Zweig-rule 
violating virtual transition i] c — > rj and r] c — > rj'. One can visualize these transitions taking 
place via the two-gluon intermediate state. Since, only the S77(3)-singlet component of the rj 
and 77' eigenstates is involved, one obtains a relation for the ratio .Rj/^ (77/77'): 



Rj/Av/v') 



T(J/jj -> r/ 7 ) 
T(J/tp -> r/7) 



1 



fc„ / tan 2 9 



(3.40) 
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where k v > and k v denote the 3-momenta of ?/ and 77, respectively. ^From the measured ratio 
(Rj/ipiv/v') = 5.0±0.8) 1461 ) a value of \9\ = 21.9° can be extracted. Thus, one gets a consistent 
result from eqs. ( |3.39| ) and ( |3.40| ). 

This Zweig-rule violating transition amplitude can also be formulated by postulating that 
the physical rj' (as well as the 77) has a small effective charm component, which should enable us 
to enact transitions of the type we are interested in, namely h — > (ccjs — > (rj, rj')s. Thus, there is 
a contribution of the operators 0\ 2 to the decay amplitude for the processes B — > (K, K*)(r)', rj), 
which can be modelled much the same way as the decays J/ip — > (t),T)')j- 

For B — > K rj', for example, this yields in the factorization approximation 

M = -^V cb V; s a 2 (v'(q)\c W c\0) (K{p')\s^b\B{p)) . (3.41) 

fc) 

The crucial quantity is the decay constant /V defined through the equation 

W(q)\ci^c\0) = -tff qfl . (3.42) 

The charm component comes in through the SU(3) singlet |r/o) which has a small charm ad- 
mixture characterized by the mixing angle 8 c5 : 

\rjo) = —j= \uu + dd + ss) cos9 c5 + \cc) sin 6^ . (3.43) 
' V3 

The orthogonal state r] c is then given by 

1 77c) = — \uu + dd + ss) sin 9 c5 + |cc) cos^g . (3.44) 

Anticipating that the mixing angle 9 c5 is small, and dropping the sin 2 9 cS term [], eq. (|3.43|) 
reads approximately as 

1 770) = —/= \uu + dd + ss) cos9 c5 + \rj c ) tan 9 c5 . (3.45) 
V3 

(c) 

/V and f Vc are then related through the equation 

ff =cos6taiL6 c£ f Tlc , (3.46) 

where f Vc is defined as (?7c(p)|c7 M 75c|0) = —if nc Pn- We estimate the r.h.s. of eq. (|3.46|) using 
experimental data. First, the mixing angle 9 c5 can be extracted from the measured ratio 

p fn/;/] J(J/^'fe7) _ (1.3±0.4)xl0- 2 ^(k v \ 3 1 

ttJ ^ nel n } ~ B( J/^ -»• 7/7) (4.31 ± 0.30) x 10- 3 " [k v ,J cos 2 ^tan 2 ^ cc ' { ' 

which on using the central values of the measurements gives \9 c5 \ = 0.014. Second, the decay 



constant f Vc can be extracted from the measured decay width f46 



4f47ra) 2 f 2 

T(Vc - 77) = l s1 ' Jvc = 7.5±H KeV , (3.48) 
8l7rm^ 



x If the mixing angle 9 c5 indeed turns out to be small, the extraction of the angle 8 discussed above is not 
significantly altered. 
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which, again for the central values, leads to f Vc = 411 MeV. Eq. (p. 48 ) is the result obtained 



using the standard nonrelativistic approach. This result also follows if one assumes that the 
shape of the distribution amplitudes for the charm and anti-charm quarks in the r] c , rf and r\ 
mesons are very similar. This gives 

\ff | = | cos 9 tan 6 c5 f Ve | ~ 5.8 MeV . (3.49) 

Similarly, we can estimate the charm content of the r\ meson: 

(V(q)\ci^ 5 c\0) = -if® g M , (3.50) 

with 

|4 C) | = |sin# tan^/J ~ 2.3 MeV . (3.51) 

Note that this method does not allow us to determine the signs of fjj, and f®, because only 
the absolute value of the mixing angle 8 c5 can be extracted. To illustrate this ambiguity in the 
numerical results, we show in the case of — > K ± r) > , the branching ratios for both signs. 

3.4.2 i] - rf — i] c system in three-angle mixing formalism 

It has been recently argued that the octet-singlet mixing scheme involving (77, 77') described 
above needs modification |4£], [5(J . More specifically, one can show that taking into account 1/N C 
corrections in the effective U(3) <8> U (3) chiral perturbation theory, one needs to distinguish the 
mixing angles of the octet and singlet components. Restricting to the (rj, rf) sector, one now 
has 

1 77) = cos6» 8 |r7 8 ) - sin^ol^o) , W) = siads\m) + cosflollo) • (3.52) 

The analogous relations for the pseudoscalar decay constants and masses can be be derived 
from the terms quadratic in = Z)o=o ... 8 -^ a 0a m the chiral lagrangian for the pseudoscalar 
mesons cj) a (here, X a ;a = 1,...,8 are the SU(3) Gell-Mann matrices and Ao is a unit matrix). 
Writing the eigenstates as <pp = J2 a /p0a> one can show that the vectors (/®, /S) and (/°, /°,), 
defined through the matrix elements involving the axial-vector currents and A Q , 

(0|4|r7(p)) = if v P„, 

miwip)} = ifipp, 

(0K\v(p)) = ifipp, 

(0|A>'(p)> = if°p^, (3.53) 



are not orthogonal to each other. Instead [49 



(fX + M) = ~^(f 2 K - fl) (1 + O(S)) , (3.54) 
where 0(8) represents 0(1/N C ) corrections. This relation then implies 

sin(# - 9 S ) = 2 ^ {Pk ~ f " ] (1 + 0(8)) , (3.55) 

J 8 
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which yields (on using the decay widths T(rj — > 27), r(77' — > 27) and the chiral perturbation 
theory result / 8 = 1.28/„-) the following values 



8 = -20.5 U , O ^ -4.0 U . (3.56) 

Thus, numerically the octet mixing angle 8$ comes out close to the angle 9 but the singlet 
mixing angle is quite small. This implies that the pseudscalar meson \t]) is almost a pure octet. 
Extending this formalism to the (77, rj', r] c ) mixing, along the lines described in the previous 
subsection, now yields the following estimates of the quantities /^f and ffl: 

flf = cos 6 tan 6 c5 f Vc 

= -sm9 tan9 c5 f Vc . (3.57) 

Using again the ratio Rj/ip{r) c /rf) given in eq. Q3.47Q yields \f$\ — 5.8 MeV (the same as 

before), but = —/if tan^o, which is considerably smaller than the previous estimate of the 
same, as 9 <C 9. 

The (77, ?/)-Hiixing framework with two angles 9q and 9 S fares better than the conventional 



one from a phenomenological point of view as well. Feldmann and Kroll [51] have compared the 



two mixing frameworks in a recent analysis of the 777 and f/7 transition form factors using data 



from CLEO J52[, L3 B TPC HI and CELLO |§, the decay widths T(rj -> 2 7 ), Ftf -> 2 7 ) 
and the ratio Rj/^(j]/r]'). They show that even after allowing for a moderate S , f/(3)-breaking 
effect, one gets a poor fit of the data in the conventional mixing formalism (i.e., with the single 
angle 9). However, the mentioned data can be well fit in the two-angle framework for the 
(77, r/)-mixing. Their best-fit values yield (with fj^ and f^P set to zero) |^TJ : 

9 8 = -22.2°, 9 = -9.1°, j- = 1.28, j- = 1.20 , (3.58) 

J IT J TV 

which agrees reasonably well with the estimates of these parameters using chiral perturbation 
theory [gj|: 

9 8 = -20.5°, 9 = -4.0°, j- = 1.28, j- = 1.25 . (3.59) 

If accurate high-g 2 data were available, one could determine the coupling constants and ffi 
from the 777 and 7/7 transition form factors, respectively. While the value fjp = is consistent 



with the data, the anaysis in |jT] yields the following range for /^f . 

- 65 MeV < < 15 MeV . (3.60) 

This determination is somewhat model-dependent as it depends on some parameters related 
to the charm wave-function. In the analysis reported in pl| , the shape of the distribution 



amplitudes corresponding to the charm quark in the r\ and 7/ are assumed to be the same as for 
r\ c . It is satisfying that the value obtained by us |/^f \ ~ 5.8 MeV from Rj/^(j] c /r]') lies within 
the range given in eq. (|3.60|) . 



In what follows, we shall adhere to the l/iV c -improved chiral perturbation theory description 
of the (77, 77')-mixing. For numerical estimates, we use the best-fit values given in eq. (|3.58|) 
above. We now discuss the decays B ± (r],r]'){K ± ,K* ± ). 
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3.4.3 B ± -> K^rf 

The matrix element M for 5~ — > reads in the factorization approximation 



.1/ 



G p 

71 



" 2 " '" i-<^ x "K) /«, 



(c) 



+ V cb V* a 2 



Jrf 

fu 
Jrf 



-VtbVZ 



+ a 4 + 



a 6 m 



2a 3 - 2a 5 + ^a 3 - a 5 + a 4 + 

2a 6 m^ \ m|-mj F W (m^) f K 



a 6 m r 



) H 

m s (m b -m s )J ft, m s (m b - m s ) 



(m s + m u ) (m b - m u ) 



m 



B 



m\ Fq~ 



■K- 



(K-\sb-\B-) (r/\uu-\0) 
(3.61) 



with 



(K-\sb.\B-) (t/|uu_|0) = ift (ml - m\) F^ K ~ (m%,) . (3.62) 

The term proportional to V cb V* s in eq. (|3.61|) is due to the charm content of the rf as discussed 
above. In eqs. ( |3.61| ) and ( |3.62| ) the decay constants ft and ft, defined as 

(0|w7 M 75«|r?') =if%Pp , (0|s7 M 7 5 s|V) = iftp^ , (3.63) 

are given in terms of f$ and fo as 



ft = -2^= sin 6» 8 + cos 6» 



(3.64) 



We remark that the matrix element (0\sj5s\r)') , which occurs when factorizing the contributions 
of O5 and Oq, has to be treated with some care. In the earlier version of the this paper we 
erroneously used the relation 

. ft m \> 



(0\sj 5 s\r)') 



2m, 



which is vitiated due to the contribution of the anomaly term in the equation 



•9^7^753 = 2m s S275S + ^G a/3 G a p 



(3.65) 



(3.66) 



To get the correct expression for the matrix element (t/|s75s|0), we now use instead the anomaly- 
free equation for the divergence of the octet axial- vector current, 

(3.67) 



(3.68) 



<9 M («7m75« + dj^d - 257^755) = 2% (m u u^u + m d drf 5 d - 2m s sj 5 i 
Neglecting the terms proportional to m u and ma on the r.h.s. of eq. ( |3.67[ ) one derives 



(Q\s~f 5 s\r)') = i 



VEf s sinflsmg, _ . (ft - ft) 7^ 



4m« 



—1 



2m, 



Of course, this relation can also be derived by working with the divergence of the (anomalous) 
singlet axial vector current. This gives rise to the term — a 6 m^/ jm s (m b — m s ) in eq. (|3.61|) . 
Likewise, the amplitudes of the other processes B~ — > K*~r]', B~ — > Kt] and B~ — > K*~r] 
published in the earlier version of this paper also get modified. The corrected amplitudes are 
given below. It appears that this (anomaly-related) error has permeated the recent literature 



73| and should be corrected accordingly. 
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3.4.4 B ± -> K^rf 

The matrix element M for — > K*~r]' reads in the factorization approximation 



M 



v/2 



F^' (m 2 K .) f. 



a 2 + ai- 



K* 



2a 3 - 2a 5 + a 3 - a 5 + a 4 - 



-v tb v* 

Ff-^Wc.) fx 



+ V cb V*a 2 



Ac) 

J y 

fu 
Jri' 



a 6 m r 



+a 4 



m s (m b + m s )J f% m s (m b + m s 
(K*-\sb_\B-) (rj'\uu-\0) 



with 



(K*-\sb-\B-) (rj'\uu_\0) = -if%2m K * (p B e^.) A^ r (mJ) 



(3.69) 
(3.70) 



3.4.5 B ± -> ^77 

The matrix element M for — > A' 77 reads in the factorization approximation 



M 



V2 



VubV* 



rnb-m* F^{ml) f K 



" r < ml Fr K -{ml) /;; 



-v tb v* 

+ I a 4 + 



a 6 m: 



f(c) 
fu 



2a 3 - 2a 5 + I a 3 - a 5 + a 4 H T/ I ^ 

\ m s (m b - m s ) J m s (m b - m s 



2a§m 2 K 



m. 



m 2 v F^\ml) f K 



(m s + m u ) (m b - m u ) J m% - m\ F^ K (777,2 ) /« 



(K-\sb-\B~) (r)\uu-\0) 
(3.71) 



where 
with 



{K-\sb4B-){ V \uu^)=i^(m 2 B -m 2 K )F^-(mi 1 ) , 



(3.72) 
(3.73) 



3.4.6 B ± -> if* ±77 

The matrix element M for B~ — > K*~r] reads in the factorization approximation 



2 + «1 



2a 3 — 2a 5 + a 3 — a 5 + a 4 



Ff^MlKl 



+ a 2 



Jri 



fs 



+a 4 



m s (m b + m s ) J f% m s (m b + m s ) 
(K*-\sb_\B-) (r]\uu-\0) 



(3.74) 
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with 



(K*-\sb„\B-)( V \uu„\Q) = -if%2m K *(p B e* K *)A%-* K (mj) . (3.75) 

4 Input Parameters, Numerical Results and Comparison 
with the CLEO Data 

4.1 Input parameters 

The matrix elements for the decay B — > h±h2 derived in the preceding section depend on the 
effective coefficients ai,...,a 6 , quark masses, various form factors, coupling constants and the 
CKM parameters. In turn, the coefficients dj and the quark masses depend on the renormaliza- 
tion scale /i and the QCD scale parameter A^jg. We have fixed A-y^r using as(Mz) = 0.118, 
which is the central value of the present world average as{Mz) = 0.118 ± 0.003 ||63|| . The scale 
p is varied between p = m h and p = rrib/2, but due to the inclusion of the NLL expressions, the 
dependence of the decay rates on p is small and hence not pursued any further. To be specific, 
we use p = 2.5 GeV in the following. The CKM matrix will be expressed in terms of the 
Wolfenstein parameters [|58|], A, A, p and the phase rj. Since the first two are well-determined 
with A = 0.81 ± 0.06, A = sm6 c = 0.2205 ± 0.0018, we fix them to their central values. The 
other two are correlated and are found to lie (at 95% C.L.) in the range 0.25 < rj < 0.52 and 
—0.25 < p < 0.35 from the CKM unitarity fits |64]]. We take four representative points in the 
allowed (p, 77) contour. Their values and the legends used in drawing the figures are as follows: 

1. p = 0.05, 77 = 0.36, yielding \J p 2 + rj 2 = 0.36 (drawn as a solid curve) 



2. p = 0.30, rj = 0.42, yielding \J p 2 + rf = 0.51 (drawn as a dashed curve) 

3. p = 0,i] = 0.22, yielding V P 1 + V 2 = 0-22 (drawn as a dashed-dotted curve) 

4. p = —0.20,77 = 0.45, yielding y/p 2 + r}' 2 = 0.49 (drawn as a dotted curve). 

They correspond to the central values of the fits in ||64j| , maximum allowed value of iKifc/Kil 
with positive p, minimum allowed value of | V^fe/V^b| , and maximum allowed value of | V^&/V^b| 
with negative p, respectively. The CKM parameters are also an output from the measured 
non-leptonic B decays and we shall illustrate the potential interest in this kind of analysis 
using some of the ratios of the branching ratios as an exercise. The rest of the input quantities 
used in our estimates for the branching ratios are collected in several tables. We discuss now 
these input values. 

4.1.1 Effective coefficients in the factorization scheme 

With the electroweak penguins and the so-called jy-annihilation/exchange diagrams neglected, 
the amplitudes for the various decays depend on six coefficients, a i; defined in section 3. Even- 
tually, one should determine each one of them (or particular combinations thereof) by analyzing 
the specific decay modes most sensitive to these coefficients. This way, one can measure the 
deviation in each one of them from their values in perturbation theory and determine if this 
deviation (due to non-perturbative effects) can be described in terms of a few universal param- 
eters. Perhaps, it should be remarked here that an analysis of the heavy to heavy transitions in 
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two-body B decays can be reasonably well described in terms of one parameter, called £ in [13], 
whose value seems to be universal. Following this, we do the simplest thing here by assuming 
that a single parameter £, defined in the preceding section, can be used to compensate for ne- 
glecting the octet-octet terms in all matrix elements of the decays B — ► h\Ji2- This is motivated 
by the fact that the energy release in these decays is comparable. It remains an open question 
if the parameter £ introduced here in the decays such as B — > Kit is close to the corresponding 
parameter £, entering, for example, in the decay B — ► Dir. We show the dependence of the 
branching ratios in several decay modes in the range < £ < 1, with £ = 1/3 being the naive 
factorization value, i.e., if one uses factorization and neglects the octet-octet contribution in 
the matrix elements. 

4.1.2 Decay coupling constants and form factors 

For the various decay constants occurring in the formulas in section 3 we use the numerical 
values shown in table 2. The values for f u , f K , j K * and f n coincide with the central values quoted 





fK 


f K * 


u 


fo 


/8 




\fV\ 


195 


160 


214 


131 


157 


168 


5.8 


0.93 



Table 2: Decay constants in MeV. 



in [|13|] extracted from data on the electromagnetic decays of u and r decays, respectively |46 
The decay constants f£, ft, /" and /* defined in eqs. ( |3.64j ) and ( |3.73| ) are obtained from the 
values for /q and fs in table 2 and by using = —22.2° for the (t],T]') mixing angle [fl7[ 
errors on the coupling constants in the table are small (typically (1 — 3)%) 



except on fy 



The 
and 



f(°\ for which present data allow a determination with an error of ±15% (assuming the mixing 
formalism holds). 

The decays being considered here, such as B — > irK, involve light hadrons in the final 
state. The rates require the knowledge of the various form factors at q 2 = m\, where 
denotes a light hadron mass. Since q 2 



mt 



is rather close to the point q = 0, and a 
simple pole model is mostly used to implement the q 2 dependence in the form factors, we shall 
neglect this g 2 -dependence in the form factors and equate F^ h (q 2 = m 2 ) = F^ h (q 2 = 0). 
Explicit calculations bear this out and find that the variation in the stated range is indeed small 



The values used for the form factors F, 



B— 
0,1 



h (q 2 = 0) and AQ^ h (q 2 = 0) in our rate 
estimates are listed in Table 3. They are taken from Q, which are reproduced in most other 
calculations (see, for example, Table 1 in fl39|l). Note also, that the 5'[/(3)-breaking effects in 
the form factors are neglected. They are typically of O(20) % 



^0,1 


TpB — >7T — 

^0,1 


TPB — >7T^ 

^0,1 


^0,1 


^0,1 


aB^uj 
^0 


aB^K* 
^0 


0.33 


0.33 


0.33 
V2 


0.33 


sin 8g 1 cos 9o 
V6 1 y/3 




0.33 


cos 6g sin do 




0.28 
V2 


0.28 



Table 3: Form factors at q 2 = 0. 
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4.1.3 Current and constituent quark masses 



The quark masses enter our analysis in two different ways. First, they occur in the amplitudes 
involving penguin loops. We treat the internal quark masses in these loops as constituent masses 
rather than current masses. For them we use the following (renormalization scale independent) 
values: 

m b = 4.88 GeV, m c = 1.5 GeV, m s = 0.5 GeV, m d = m u = 0.2 GeV. (4.1) 

Variation in a reasonable range of these parameters does not change the numerical results of 
the branching ratios significantly. The value of nib above is fixed to be the current quark mass 
value Tribal = m&/2) = 4.88 GeV, given below. Second, the quark masses nib, rn s , ^d and 
m u also appear through the equations of motion when working out the (factorized) hadronic 
matrix elements. In this case, the quark masses should be interpreted as current masses. Using 



nibinib) = 4.45 GeV 66 and 



mj(l GeV) = 150 MeV , rfQ{l GeV) = 9.3 MeV , m^(l GeV) = 5.1 MeV , (4.2) 



from ||62|| , the corresponding values at the renormalization scale \i = 2.5 GeV are given in table 



4, together with other input parameters needed for our analysis. 



m h 


m s 


m d 


m u 


a s (m z ) 


TB 


(d 2 ) 


08 


0o 


4.88 GeV 


122 MeV 


7.6 MeV 


4.2 MeV 


0.118 


1.60 ps 


mi/2 


-22.2° 


-9.1° 



Table 4: Quark masses and other input parameters. The running masses are given at the 
renormalization scale fi = 2.5 GeV. 



4.1.4 Numerical values for the effective Wilson coefficients 

From eqs. ( |2.jj| ) - ( ^.10[ ) it follows that the effective Wilson coefficients defined in eq. ( |2.5| ) 
are in general complex numbers, which depend on quarks masses and on the CKM matrix 
elements. Taking the quark masses listed in eq. ( [4.1|) and using the central values for the 
CKM parameters from the unitarity fits |64| (i.e. p = 0.05 and t] = 0.36), the effective Wilson 
coefficients Cf at the renormalization scale p = 2.5 GeV are shown in table 5. We remark 
that the (almost) identical values of these coefficients in the first two columns (b — > s and 
b — > s) reflects that the imaginary parts of these effective Wilson coefficients are essentially 
generated by strong interactions. The numerically differing entries in the other two columns 
(b — > d and b — > d) reflect that the weak (CP-violating) and strong interaction phases in these 
decays are comparable. 

4.2 Numerical results and comparison with CLEO data 

Having stated our theoretical framework and the input parameters, we now present our results 
for the various decays of interest listed in the previous section. A word of caution concerning 
the accuracy of the absolute decay rates calculated by us is in order. As just displayed, there 
are many parameters involved in describing exclusive non-leptonic decays and while the decay 
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h — > rl 


h — > rl 


(jet 3 


1.160 


1.160 


1.160 


1.160 


c%" 


-0.334 


-0.334 


-0.334 


-0.334 


c e / f 


0.021 + 0.004i 


0.021 + 0.004i 


0.020 + 0.002i 


0.022 + 0.006i 


c\" 


-0.052 - O.Olli 


-0.051 - O.Olli 


-0.048 - 0.007i 


-0.053 - 0.017« 


cf f 


0.016 + 0.004i 


0.016 + 0.004i 


0.015 + 0.002i 


0.017 + 0.006i 


ci" 


-0.064 -O.Olli 


-0.063 -O.Olli 


-0.060 - 0.007i 


-0.065 - 0.017z 



Table 5: Effective Wilson coefficients C\ at the renormalization scale \i = 2.5 GeV for the 
various b — > q (b — > q) transitions. See text and eq. Q2.8Q . 

rates do not depend sensitively on all of them, and many input parameters are already well 
known, it is obvious that the predicted branching ratios do depend sensitively on some for 
which there is no alternative at present to using model-dependent estimates. The particular 
quantities in question are the decay form factors. Some of these form factors enter in other 
processes which have been measured (such as in the semileptonic and radiative B decays) and 
the estimates being used are found to reproduce the data quite well, yet some others are not yet 
constrained by data directly. So, the estimates given below for the absolute decay rates have to 
be taken with an accuracy which is not better than a factor 2. The additional uncertainty due 
to the parameter £ can not be judged at this stage. That can only be ascertained in future, if 
this framework proves to be a reasonable way to analyze heavy to light transitions in B decays. 

However, within this framework, the ratios of the branching ratios are much more stable, 
as many of the theoretical uncertainties (such as in the form factors, various scales, and quark 
masses) cancel out to a large extent. In some cases, the dependence on the parameter £ also 
cancels, or it is very weak. Hence, the ratios are more reliable and the experimental information 
on these ratios can eventually be used meaningfully to draw inferences on the fundamental 
parameters, such as p and rj. 

4.2.1 Branching ratios for B — > tttt modes 

We shall show the branching ratios of interest as a function of the parameter £ for four different 
set of values of the CKM parameters. Wherever available, the present measurements of the 
branching ratios at the ±1<t level are also shown on these figures (thick solid lines). All ex- 
perimental numbers are taken from [|Ij, El H], and in showing the experimental results, we have 
added the statistical and systematic errors in quadrature. We start by showing in Fig. ^ the 
branching ratio B(B ± — > 7r°7r :t ). The decay rate for this mode is sensitive to both the variation 
in £ and the CKM parameters. This is obvious from the quadratic dependence of the decay 
rate on the quantity \V u b\. Also, it depends on the combination ci\ + ci2- Hence, a measurement 
of this decay rate will yield information on these quantities. In quoting a range, we shall take 
< £ < 0.5 (which is suggested by the combined analysis of all the present CLEO data on 
B — > hxhz decays which we show later). We estimate, 

B{B ± -> ttV) ~ (0.1 - 1.4) x 10~ 5 , 
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which is uncertain by over an order of magnitude. However, the lower range corresponds to 
the rather small value of the CKM-factor, IV^/V^I = 0.05, and is therefore somewhat unlikely. 
For the central value \V u b/V c b\ = 0.08, we estimate 

B{B ± -> vr ^) = (0.3 - 0.6) x 10~ 5 . 

The present experimental upper limit is (at 90% C.L.) 

B{B ± -> ttV) < 2.0 x 10~ 5 . 

In Fig. |5], we show the branching ratio B(B°(B°) — >• 7r =F 7r ± ). Again, this decay mode is 
sensitive to £ and the CKM parameters, although the resulting uncertainty is less in this case 
than in B(B ± — > 7r°7r ± ). Comparison of the model calculations with the present upper limit 
(at 90% C.L.) 

B(B°(B°) -> ttV) < 1.5 x 10" 5 . 

shows that this decay mode is expected to lie within a factor 2 — 3 of the present upper limit 
and hence should be measured soon. Already, the present upper limit on this mode disfavors 
some extreme values of the CKM parameters corresponding to \ V u b/V c b\ close to or in excess of 

0.11. 

In Fig. |6], we show the branching ratio B(B°(B°) — > 7r°7r°). This branching ratio is not very 
sensitive to £ in the region < £ < 0.5 but rises sharply as £ — ► 1. All the curves lie however 
significantly below the present upper limit p| : 



B(B°(B°) -» ttV) < 9.1 x 10~ 6 . 
Restricting to < £ < 0.5, our model calculation yields 

B(B°(W) -> ttV) ~ (0.5 - 2.0) x 10~ 6 . 

4.2.2 Branching ratios for B — ► ixK modes 

In Fig. [7|, we show the branching ratio B(B ± — > 7T K). This is a good decay mode, in principle, 
to determine the parameter £, as there is no perceptible dependence of the rate on the CKM 
parameters. In the indicated range < £ < 1, the branching ratio varies by slightly more than 
a factor 2. The experimental measurement is (at ±3.2<r): 

B(B ± -> ^K) = (2.3±i:J±g;l ± 0.2) x lO" 5 . 

Our estimated branching ratio is in agreement with data, and there is a slight preference for 
smaller values of £, with £ > 0.7 somewhat disfavored. Since the CKM-parametric dependence 
is small, this decay mode is useful to show the effects of the QCD corrections. In Fig. [5], we show 
the branching ratio B(B ± —>■ ^K) as a function of the scaled variable (q 2 /ml), in the range 
< (q 2 /iTil) < 1, calculated for £ = 0. The dashed line corresponds to the LL approximation, 
whereas the dotted and solid lines correspond to the truncated NLL approximation, and the 
complete NLL approximation as discussed in section 2, respectively. The dotted curve amounts 
to what has been used in the analysis of the decay modes (B — > tt^K) in |2~8 , f2~9fl . The 
effect of the complete NLL corrections is numerically important, and they tend to decrease the 
branching ratio as compared to what one estimates by including the charm penguins alone. 
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In Fig. we show the branching ratio B(B°(B°) — > ti ± K t ). Like its charged partner, 
B(B ± — > -n^K) discussed above, this decay mode is also sensitive to the parameter £, though 
in this case there is a perceptible dependence of the rate on the CKM parameters as well. The 
observed branching ratio (at ±5.6<r): 

B(B°(W) -> ttKK^) = (1.5+0 5+0.1 ± o X ) x 10 -5 9 

is quite comfortably accommodated by our estimates. 

Comparing Figs. |9| and |7|, one sees that the dependence of these decay rates on £ is very 
similar, and hence in the ratio of branching ratios it almost cancels out. Defining this ratio by 
R 

_ B(B°m -> «±K*) 
Rl = B(B±^n±K) ' (43) 

we show Ri as a function of the CKM parameter p in Fig. [IJJ for two values of the CKM 
parameter rj = 0.52 (upper curve) and rj = 0.25 (lower curve). We note that Ri is rather 
insensitive to r\ but it does depend sensitively on p. Using the present CLEO measurement of 
Ri (at ±la) 

Ri = 0.65 ± 0.40 , (4.4) 



Fig. [10] suggests that negative values of p are disfavored. This can also be converted as a 
statement on the CP-violating phase 7. Since the Wolfenstein parameter 7] is positive as 
determined from the constraint on ex, p > implies 7 < 90°. We recall that the bounds on 
7 obtained from the CKM unitarity fits yield symmetric constraints, centred around 7 = 90° 
(or p = 0). However, it should be remarked that the lower bound on the ratio of the weak 
mass differences in the B° s -B® and B®-B® systems, AM S / AM d , which at present is posted as 
AM S / AM d > 20.4 at 95% C.L. p7| , now cuts away a good part of the negative- p (equivalently 



7 > 90°) region. A recent analysis gives (at 95% C.L.): 32° < 7 < 122° |68| , which is no 
longer symmetric around 7 = 90°. On the other hand, the model-independent constraints on 
7 from Ri, discussed by Fleischer and Mannel p8 |, are such that they force 7 to lie in the 



range 0° < 7 < ^ max or 180° — j max < 7 < 180°, depending on the sign of cos 5, where 
5 is the strong phase-shift difference between the tree and penguin amplitudes in the decay 
B°(B°) — > tt ± K t . Since this phase difference is calculated in our model, the preferred solution 
is the one in which 7 lies in the first quadrant, or 0° < 7 < ^ max . Unfortunately, with the 
present experimental errors, the 95% C.L. limit on p from Ri (or on ^ max ^ does not allow one 
to draw more quantitative conclusions on the value of 7 than what one gets from the CKM fits 
pi, Blfl . This is expected to change with improved data on Ri, if the value of Ri is found to be 



considerably less than 1. Our analysis, carried out in the factorization framework, underlines 
the sensitive dependence of Ri on p, with p < disfavored (at ±1<t) by the CLEO data on R\. 
The effect of the present lower bound on AM S / AM d on p is qualitatively similar to the one 
from the present measurement of R\, namely both prefer p > 0. For an updated CKM fits, see 



also 65 . 



4.2.3 Branching ratios for the B — > fv^n and B — > h^K modes 

The decay modes B ± — > h ±r K° , B°(B°) — > h ±r K T and B ± — > h^K have been measured with 
impressive precisions. We compare our model estimates with these measurements. In Fig. |lj], 
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we show the branching ratio B(B ± — > h 7T°). The decay rate in this case is mildly dependent on 
£, but more importantly on the CKM parameters. The experimental measurement (at ±5.5ct): 

B(B ± -> tt ^) = (1.6±g;i±gJ ± 0.1) x 1(T 5 

is reproduced well by our model. 

In Fig. [T2], we compare our model estimates with the CLEO measurements (at ±7.8ct): 

B(B°{W) -> h***) = (2.2t°j ± 0.1) x 10~ 5 . 

Agreement between our model and data is good. The two curves (dashed and dotted) which lie 
outside the ±lcr bands correspond to large values of the ratio \V u b/V c b\, namely |Kife/Kb| = 0.11, 
which is also outside of the ±1ct bound from direct measurements of IV^/V^I- So, all of these 
different pieces of data are giving a consistent picture. 

In Fig. 13, we show our estimates for the branching ratio for the mode B ± — > h^K , which 
has been measured (at 4.4cx) 

B(B ± -> ^K ) = (2.4±i;5±g;l ± 0.2) x 10~ 5 . 

This branching ratio has a very similar dependence on £ as in the decay B ± — > n^K and 
likewise has little dependence on the CKM parameters. Model estimates are in agreement with 
data for £ < 0.7. 

As another example of a ratio of branching ratios, which is sensitive to the CKM parameters, 
we define the ratio R2 

_ B (B°m -> g7g) 

which like Ri is less dependent on the other input parameters, including £. Since B(B ± — > 
^K ) is insensitive to the CKM parameters, the ratio R2 reflects the CKM dependence of 
B(B°(B°) — >• /i ± 7r T ). We plot the ratio R2 in Fig. [TJ| as a function of the phase 7], for three 
values of p: p = 0.05 (dashed curve), p = 0.35 (solid curve) and p = —0.25, which coincides 
with the case p = 0.35. The present experimental value of R2 (at ±la) is 

R 2 = 0.96 ± 0.57 . (4.6) 

This shows that with the stated significance R2 disfavors large values of t] in excess of rj > 0.5. 



4.2.4 Branching ratios for the B ± — >• ujK ± and B^ — >• tuh^ 1 modes 

Next, we study the decays B ± — > ujK ± and B^ — > ioh^ (h = n,K), which have also been 
measured by the CLEO collaboration with the former having a branching ratios (at ±3.3<r) 

B(B ± -> uK*) = (1.2^ ± 0.2) x 10~ 5 , 

and the latter (at ±6. Oct) 

B(B ± -> uh*) = (2.5t° 8 7 ± 0.5) x 10~ 5 . 

These measurements are compared with our model calculations in Figs. [15] and [16], respectively. 
Both of these decays have an interesting dependence on the variable £. Taken the data at face 
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value (±lcr), a value for £ in the range 0.15 < £ < 0.5 and £ > 0.85 are somewhat disfavored 
by data in the decay B ± — > ouK ± . Curiously, the estimated branching ratio B{B ± — > ujK^) 
has its lowest value in the range £ = 0.3 ± 0.1, and in this range it fails to reproduce the data 
by almost 2a. This observation and the present measurement of B(B ± — ► tt^K) as well as 
B{B ± —>■ h^K), which disfavor £ > 0.7 then imply that the preferred value of £ in our model 
is either in the range < £ < 0.15, or else £ ~ 0.5. In this range, however, the estimated 
branching ratio is somewhat lower than the experimental one in B ± — ► uih , but not by a large 
amount. Due to the fact that the data being discussed are the first ones of their kind and 
the uncertainties related to the parameters of the present theoretical framework have not been 
exhaustively studied, one can not draw too strong conclusions on the value of the parameter £ 
from this decay. 



4.2.5 Branching ratios for the B ± — > (rj,r)')(K ± , K^) modes 

Finally, we take up the decay B — > rj'K , which has attracted a lot of theoretical attention 
recently. Compared to the decays considered so far, this decay and the related ones B ± — > r]K ± , 
B ± — ► ri'K*^ 1 and B — > r)K have an extra contribution from the decay chain b — > sec — ► 



In Fig. [17] we show the branching ratio B{B — > q'K ) as a function of £, varying the 



CKM parameters as indicated in section 4.1. Since we are not able to determine the sign of the 
coupling constant fty due to the sign ambiguity in the determination of the angle 6 C , we show 

(c) 

the result for both f^, = ±5.8 MeV. Note that the ^-dependence of this branching ratio results 
in a factor 2 uncertainty varying £ in the range < £ < 0.5 for the positive-/^/ solution; the 

branching ratio is less sensitive to £ for the negative-/^, case. The positive-/,^ solution yields 
a marginally higher branching ratio. The CKM-parametric dependence of this branching ratio 
is not very marked. Within the present uncertainties in the input parameters, we get, at £ ~ 

B(B ± -> ri'K*) ~ (3 - 4) x 10~ 5 , 
which at £ = 0.5 falls down to the range 

B{B ± -> n'K*) ~ (2 - 3) x 10" 5 . 

This is to be compared with the CLEO measurement (at ±5. 5c) 

B{B ± -> = (7.1±l;f ± 0.9) x 10" 5 . 

Given the experimental and theoretical errors, the model estimates and data are clearly not 
incompatible though, with the values of the parameters used by us, our estimates are somewhat 
on the lower side. Since, apart from the form factors, this branching ratio is sensitive to the 
value of the s-quark mass, with the branching ratio increasing as m s is decreased from its 
default value m s (fj, = 2.5 GeV) = 122 MeV used by us, the decay rate can be made to agree 
with the CLEO data by optimising these parameters in an overall fit. This is not warranted at 
this stage. 

The branching ratios for the decays B ± — rj'K*^, B ± —> r]K ± and B ± — >• r)K are shown 
in Figs. 0, 0, and |0], respectively, for the values f { f = -5.8 MeV and = -0.93 MeV. The 



reason for selecting the negative sign is to be traced to the observation that the contribution 
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of the amplitude b — > (cc)s — > (j],r]')s can also be calculated using QCD-anomaly which fixes 



the signs of these constants to be negative |69]. The estimated branching ratios satisfy the 



respective present experimental bounds on them [|IJ. For < £ < 0.5, we predict 

B(B ± -> ri'K**) ~ (0.3 - 0.9) x 10~ 6 . 

The decays .B^ — > ri(K ± ,K ± *) on the other hand do not depend very sensitively on the sign 
of /W. We estimate (0 < £ < 0.5) 

B(B ± -> 77^) ~ (1.0 - 2.8) x 10~ 6 , 

-> rjK^) ~ (1.0 - 2.8) x 10~ 6 . 

Finally, we remark that scenarios with a greatly enhanced strength of the dipole operator 
Og have been entertained in the literature j|9], |6(| plR , with the view of bringing the existing 
theoretical estimates of the semileptonic branching ratios and charm counting in B decays in 
better rapport with data. A greatly enhanced value of C&(mw) will influence the branching 
ratios in some selected non-leptonic B decay channel as well. However, this effect is diluted due 
to the contributions from other Wilson coefficients, which are assumed to have their SM values. 
Also, as emphasized in [[5^, |5D[ , and more recently in |plf , the strong mixing of the operators 
O2 and 0$ would require a very large enhancement in C%(mw) I 'Cf M (mw) , typically O(10), to 
have a measurable influence in B decays, calculated at the scale \i ~ m&, due to the effects of 
the renormalization group. Qualitatively, this picture also holds in the analysis of the exclusive 
non-leptonic B decays discussed by us. We show a typical case B ± -> Kir* in Fig. |l|, where 
the branching ratio for this mode is plotted as a function of the variable C 8 (mw) / 'Cf M (mw) ■ 
Despite the large range of this variable, we find that the influence of such a markedly enhanced 
Cs(mw) on non-leptonic B decays is marginal. In future, we hope that these matters will be 
scrutinized much more minutely. We conclude that the non-leptonic -B-decays considered here 
do not require large enhancements of Cg(mjy), or of any other Wilson coefficient, as they are 
by and large compatible with data with their SM values. 



5 Summary 

In the first part of this paper we have presented a theoretical framework to study two-body 
decays of B mesons with two light mesons in the final state. First, we took into account the 
complete NLL corrections at the partonic level, improving thereby previous calculations. In 
particular, we have also included the effects of the chromomagnetic penguin operator in non- 
leptonic B decays. These NLL order corrections are numerically important in the exclusive 
decay rates. Second, to estimate the hadronic matrix elements we assumed factorization and 
gave a parametrization for the so-called heavy-to-light transitions. In its most economic version, 
this brings in a single phenomenological parameter, called £, which has to be determined by 
comparing the predictions of this model with data. We have taken a first look at the available 
CLEO data and estimated that our model is compatible within the present theoretical and 
experimental errors with data in the range < £ < 0.5. With more precise data one should 
be able to test our model and see if within reasonable accuracy one obtains a universal value 
for this quantity in the heavy-to-light sector. Alternatively, with more precise data in several 
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decay modes, we propose to extract the effective coefficients a\, a 6 directly to determine the 
extent of non-perturbative effects in each one of them. 

In the second part of this paper we have applied this framework in the analysis of the exclu- 
sive two-body B decays, in which QCD penguins play an important role. Some of these decays 
have been measured recently by the CLEO collaboration 0, [J , with which we compared our 
model calculations; we have also predicted the branching ratios for some related decay modes 
which have not been measured yet. While the formalism provided here is generally applicable 
to study all B — ► PP, B — > PV and B — > VV decays, we have restricted ourselves to discussing 
the four generic cases: B — > tttt, B — > Kw, B ± — > Ldh^ 1 and B ± — > (rj, r] > )(K ± , K*^). 

In particular, we have studied at some length the last class, involving the decay B — > K rj 
and the related ones. As the rj and r\ mesons are composed of u, d and s quarks, the correspond- 
ing decay rates are particularly sensitive to interference effects among the several competing 
amplitudes involving the current- current and the QCD-penguins operators, as was pointed out 
by Lipkin some time ago ||. In addition, the operators 0\ 2 which induce transitions of the 
form b — > s(cc) — > s(i],i]') have to be included. Estimates of the latter require a trustworthy 
evaluation of the cc component in the wave-function of the rj and rj mesons. We have used 
the mixing formalism involving the (rj, rj', ?] c )-complex and data to determine the cc contents of 
these mesons. We find that this charm-induced contribution does not dominate the amplitudes 
for the processes involving rj\ the decay rate is more sensitive to the penguin contributions. 
Our paper provides the complete amplitudes showing all these individual contributions. This 
can be used in future analyses of more precise data to determine the cc components in r\ and rj . 
The estimates presented here with |/£ c *| = 6 MeV yield B{B ± -> K±rj) = (2 - 4) x 10~ 5 . This 
is somewhat lower than the central value of the present measurement but compatible with the 
value obtained by fluctuating down the experimental error by la. However, a simple answer 
about the large measured value of B(B ± — > rjK ± ), a question frequently asked, in terms of 
a single dominating amplitude is not readily available, though the penguins and the singlet 
component of the rj are certainly at the back of the enhanced branching ratio for this decay. In 
our analysis, we find that the measured rate in the r)'K ± mode is only marginally (say a factor 
2) larger than our model estimates and given the theoretical errors there is nothing anoma- 
lous about it. We also expect that the data will evolve with time so as to reduce the present 
discrepancy. 

We have made predictions to test this interference pattern in the related decays involving i] 
and rj . The resulting decay rates, which also reflect the built-in angular momentum differences 
between the states K*(j], rj') and K(j], r]'), show a certain hierarchy among the branching ratios. 
While the other three may turn out comparable with each other (within a factor 2 -3), we 
predict: 

B{B ± -> rj'K*) > i3(fi ± -> ri'K**). 

The measurement of B{B ± — > t]'K ± ) being the largest measured so far is in line with our analysis 
based on the SM. However, in our SM-based framework it would be difficult to accommodate 
a much larger branching ratio B(B ± — > rjK ) vitiating this hierarchy. 

The rates for the other decays presented in our analysis are also in reasonable agreement with 
data, within the presently allowed CKM-parameter space. Based on our estimates presented 
here, we expect the decay mode B°(B°) — > n T n ± to be measured within a factor 2-3 below 
the present upper limit. We point out interesting inferences which present data allows to draw 
on the consistency of the SM. In particular, the ratios R\ and R2 involving the Kir and tttt final 
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states appear very promising. Present measurements on these ratios are tantalizingly close to 
providing independent information on the CKM-Wolfenstein parameters p and 77. Their impact 
on the CKM phenomenology will only be determined with more precise data, to which we look 
forward with animated interest. For the time being, the standard model rules OK - also in the 
non-leptonic B decays! 
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Figure 4: Branching ratio for B — > tt ^ as a function of £ for various points in the (p; 77)- 
plane. Solid curve: (p, 77) = (0.05,0.36); dotted curve: (p, 77) = (—0.20,0.45); dashed curve: 
(p, 77) = (0.30,0.42); dashed-dotted curve: (p, 77) = (0.00,0.22). The horizontal thick solid line 
(with the arrow) shows the CLEO upper limit (at 90% C.L.). 




34 



O 



0.8 



n — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 




Figure 6: As Fig. |] but for the process B°(B°) — > n n° 
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Figure 7: Branching ratio for — > ir^K as a function of £ for various points in the (p, 77)- 
plane. Solid line: (p, rj) = (0.05,0.36); dotted line: (p, 77) = (—0.20,0.45); dashed line: (yo, 77) = 
(0.30,0.42); dash-dotted line: (p,r?) = (0.00,0.22). The thick solid lines show the CLEO 
measurement (with ±lcr errors). 
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Figure 8: B(B^ — > ir^K) as a function of (q 2 ) /ml. The dashed line corresponds to the LL 
approximation. The solid line and the dotted line correspond both to the Wilson coefficients 
evaluated in the NLL approximation; the solid line takes into account the penguin diagrams 
of all the four- Fermi operators and the tree level matrix element of Og, while the dotted line 
takes into account the penguin diagrams associated with the four- Fermi operators 0\ 2 only. 
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Figure 9: As Fig. |7| but for the process B°(B°) — > ir^K^. 



36 



2 



1.5 

K 1 

0.5 





_1 I I I I I L_ 



0.2 



0.2 
P 



Figure 10: The ratio R x = B(B°(B°) -> n ± K T ) /B(B ± -> ifv^) as a function of the Wolfen- 
stein parameter p, for rj = 0.25 (lower curve) and r] = 0.52 (upper curve). The curves are drawn 
for £ = 0. The horizontal thick solid lines show the CLEO measurement (with ±lcr errors). 




Figure 11: As Fig. [7|but for the process B ± — * h^ir (h = K, 
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Figure 14: The ratio R 2 = B{B°{B°) -> h ± 7i^)/B{B ± -> /frr*) as a function of the Wolfen- 
stein parameter i], for p = 0.05 (dashed curve) and p = 0.35 (solid curve). For p = —0.25 
the corresponding curve is almost identical to the curve for p = 0.35 The curves are drawn for 
£ = 0. The horizontal thick solid lines show the CLEO measurement (with ±lcx errors). 




Figure 15: As Fig. [7| but for the process B ± — > u>K ± . 
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Figure 16: As Fig. [7] but for the process B ± — > cu/i ± . 




Figure 17: Branching ratio for B ± — > r]'K ± as a function of £ for the same points in the (p, rf)- 
plane as in Fig. 0. The upper (lower) set of curves close to £ = corresponds to the positive 
(negative) solution for fjj\ The horizontal thick solid lines show the CLEO measurement (with 
±lcr errors). 
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Figure 18: Branching ratio for B ± — > rj'K^* as a function of £ for the same points in the 
G°> ^)-plane as in Fig. [7]. All curves correspond to the value f$ = —5.8 MeV. The upper limit 
from CLEO is 2.9 x 10" 4 at 90% C.L.. 
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Figure 19: Branching ratio for B — > 7]K ± as a function of £ for the same points in the 

0- 

from CLEO is 8 x 10~ 6 at 90% C.L.. 



(p, ?7)-plane as in Fig. [j]. All curves correspond to the value fff' = —0.93 MeV. The upper limit 
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Figure 20: Branching ratio for B ± — > rjK^* as a function of £ for the same points in the 
(p, r^-plane as in Fig. [7|. All curves correspond to the value f( c > = —0.93 MeV. The upper limit 
from CLEO is 2.4 x 10" 4 at 90% C.L.. 
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Figure 21: B(B ± — > K^) as a function of C%(mw) I 'Cf M (mw) , where C$(mw) is the Wilson 
coefficient of 0$ including new physics, while C$ M (mw) is the standard model value. We use 
f = and (p,rj) = (0.05,0.36). 
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